Plate Boundary Research Articles

Neotectonics on the Namuncurá (Burdwood) Bank: Unveiling seafloor strike-slip processes along the North Scotia Ridge

The North Scotia Ridge is the offshore morphostructural expression of the left-lateral transcurrent South America-Scotia Plate Boundary. Several blocks compose the ridge, including the scarcely studied Namuncurá Bank (NB, also known as Burdwood). We present the first detailed study of active structures on the seafloor of the western NB from a database of 3D and 2D seismic data, multibeam bathymetry, and sub-bottom profiles. This work assesses the architecture, style of deformation and Cenozoic evolution of NB, where several groups of faults, and en echelon folding, affects the seabed and shallow sub-bottom. These features compound the northernmost structures associated with a releasing bend, fitting well with a left-lateral Riedel shear model oriented at N74°E, slightly rotated with respect to the present-day plate boundary stress regime. The current tectonic scenario started with a main deformational phase in the Neogene, partially distributed by the Malvinas Fold-Thrust Belt, while modern deformation continues to be conditioned by pre-existing structures. This study allows for a better understanding of the tectonics of the North Scotia Ridge, a morphostructure that influences the circulation of the Antarctic Circumpolar Current thus impacting the global climate.

Along‐strike variation of fault slip rate of a transform plate boundary in Tierra del Fuego (South Patagonia)

AbstractThe Magallanes–Fagnano Fault is an active strike‐slip structure accommodating the relative displacement between South America and the Scotia Plate. The Chilean portion of the fault is poorly studied because most of it runs below the sea level in the Strait of Magellan. Our tectonic geomorphological study is focused on a rare onshore fault section, along which streams horizontally deflected by hundreds of metres since the last main deglaciation are compatible with a dominant left‐lateral fault kinematics and yield a slip rate of 15.7 ± 2.4 mm/year. This rate is between 2 and 3 times higher than the one estimated on the Argentinian portion of the fault over the same period. This spatial variation may be due to both glacial unloading on the fault zone and/or structural factors. These results point out the need to study strike‐slip faults on several portions to unravel behaviour changes related to internal or external forcing.

Double-Vergent Plate Boundary Faults and Triggered Coseismic Rupture of the 2022 Chihshang Doublet Earthquake Occurred in Eastern Taiwan

Abstract The earthquake doublet of Mw 6.5 and 6.9 occurred along the west-dipping Central Range fault (CRF) adjacent to the east-dipping creep segment of the Longitudinal Valley fault in eastern Taiwan on 17–18 September 2022. The faulting model derived from the Interferometric Synthetic Aperture Radar and Global Positioning System observation suggests that the west-dipping CRF rupture is responsible for the two mainshocks. Meanwhile, the two major earthquakes resulted in ∼100 km of fault slip along the double-convergent plate boundary faults near the Longitudinal Valley in eastern Taiwan. The 2022 Chihshang earthquake sequence filled the seismic gap of the CRF located between the 2006 ML 6.1 Taitung earthquake and the 2013 ML 6.1 Ruisui earthquake. Finally, the 2022 Chihshang earthquake sequence increased the Coulomb failure stress in the southernmost segment of the CRF, which may cause a higher seismic potential in the future.

Unfolding rotational tectonics and topographic evolution from localized verses diffuse plate boundary counterparts

We present a kinematic model developed from geodetic observations, topography analysis and analogue tectonic modelling results, which reveals a striking similarity between the rotational tectonic settings of the Gakkel Ridge-Chersky Range system in the Arctic, and the Central Indian Tectonic Zone within the Indian subcontinent. A crucial aspect of large-scale extensional rift systems is the gradual variation of extension along the rift axis, due to plate rotation about a Euler pole, which may lead to contraction on the opposite side of the Euler pole to form a rotational tectonic system. Our geodetic and topographic analysis, combined with the reanalysis of analogue tectonic modelling results demonstrates such rotational tectonic plate motion in both the Arctic and Indian case. However, the plate boundary between the North American and Eurasian Plates as represented by the Arctic Gakkel Ridge-Chersky Range system is strongly localized, whereas the Central Indian Tectonic Zone that separates the North and South India Plates involves diffuse deformation instead. Furthermore, in both the Arctic and Central Indian we find that the relative Euler rotation pole is located near an indenter-like feature, which possibly controls the present-day rotational tectonics and contrasting topography on opposite sides of the Euler pole.

Recent advances in characterizing the crustal stress field and future applications of stress data: Perspectives from North America

Abstract The stress field controls patterns of crustal deformation, including which faults are likeliest to cause earthquakes or transmit fluids. Since the 1950s, maps of maximum horizontal stress ( S Hmax ) orientations have advanced dramatically, and the style of faulting (relative principal stress magnitudes) has recently been mapped in some regions as well. This perspectives paper summarizes developments in characterizing stress orientations and (relative) magnitudes, including new seismic and borehole methods, as well as progress in identifying the causes of stress variations. Despite these advances, adding far more spatiotemporal detail would allow geoscientists to address many of today's key challenges regarding natural hazards, energy development, and geodynamics. In particular, it is critically important to characterize stress heterogeneity at multiple scales while also recognizing the coherent variability of the stress field. The second part of the paper considers how more detailed stress datasets could prove essential to addressing some of the grand questions in geoscience, including deciphering the poorly understood feedbacks between crustal dynamics and surface processes, improving earthquake and eruption forecasts, and determining the origins and shared properties of plate boundaries.

Dynamical Modeling of Fault Slip Rates at the New Zealand Plate Boundary Indicates Fault Weakness

AbstractWe construct a thin‐sheet dynamical model of the New Zealand plate boundary that includes faults. Our model fits fault slip rates, style of distributed deformation, and is constrained by relative plate boundary motion. We assume a pseudo‐plastic rheology and achieve a best fit to slip rate observations with a deviatoric stress magnitude of 20 MPa. Modeled local forces are significant at Puysegur and Hikurangi subduction zones, and smaller forces are related to mantle downwelling beneath South Island and Havre Trough mantle upwelling. Modeled tractions on faults are mostly 5–20 MPa, similar to or slightly smaller than stress magnitudes adjacent to faults. Modeled shear tractions are generally 2–10 MPa, comparable to stress drops during earthquakes. Modeled stress orientations and fault dips suggest that many faults are not optimally oriented for their style of faulting. Notably small traction magnitudes of <5 MPa and shear tractions of <0.5 MPa are modeled for faults in the central North Island Dextral Fault Belt (NIDFB), which we infer to be very poorly oriented. Friction coefficients on faults (ratio of shear stress to effective normal stress) are in the range 0.1–0.3 for major crustal faults such as the Alpine Fault and Marlborough faults, but subduction zones and the NIDFB have values <0.1. We propose that low values of long‐term fault strength, shear stress resolved onto the fault, and overall magnitudes of deviatoric stress in the crust are a consequence of dynamic weakening of faults during fault slip.

Robust sulfur and oxygen isotope evidence for a highly anoxic paleoenvironment in Late Devonian seawater: Insights from marine anhydrites in the Zaige Formation, South China

A global-scale mass extinction event occurred in the Late Devonian, and its triggering mechanisms have remained a subject of controversy. The key to resolving this controversy lies in elucidating paleoenvironmental characteristics and their coupling relationship with the mass extinction. Marine evaporites serve as sensitive indicators that can be used for reconstructing the geochemical history of ancient oceans. In this study, we present a dual sulfur and oxygen isotopic dataset of marine anhydrites in the Lunatian gypsum deposit (Yunnan, China), and investigate its paleoenvironmental significance within the framework of the Late Devonian mass extinction. These anhydrites originate in arid and hot continental margin seas near convergent plate boundaries, forming as a result of the dehydration of marine gypsums during burial processes. Anhydrites exhibit a conspicuous positive anomaly in δ34S values (23.6–25.1‰, averaging 24.3‰), primarily attributed to bacteria-mediated sulfate reduction. The δ18O values of the Lunatian anhydrites (12.2–15.2‰, averaging 12.0‰) indicate a slightly higher oxygen isotope composition for Late Devonian seawater (8.7–11.7‰, averaging 10.5‰) compared to those of contemporary seawater (9.5–10.1‰). Such an oxygen isotopic characteristic for the Lunatian anhydrites is intricately associated with the chemical processes that occur during the diffusion of sulfides into the upper seawater. The anhydrite-forming evaporating system has a noticeable brine layering where the surface seawater is less dense than the deeper brine, and bacterial processes gradually increase δ34S in the deep brine. When the deep brine and surface seawater rapidly mix, it leads to a significant positive shift in the sulfur isotope composition of the resulting evaporites. Our sulfur and oxygen isotope data from Lunatian anhydrites validate the existence of an exceedingly anoxic marine environment during the Late Devonian, potentially serving as a pivotal factor in the mass extinction event.

Continued Convergence After the Occurrence of a Slab Break‐Off: The Case of the Cyprian Arc

AbstractThe detachment (i.e., break‐off) of down‐going subducting oceanic slabs is a major geodynamic event with far‐reaching consequences, one of which is the reduction of the slab pull force acting on the trailing plate. We investigate the motion of the Sinai Microplate where a recent (∼1 Myr ago) slab break‐off occurred along its sole converging plate boundary (Cyprian Arc) with the overriding Anatolia Microplate. Based on new bathymetric mapping, high‐resolution seismic reflection imaging, geodetic and earthquake data, we show that Sinai is actively moving in a northwest direction with respect to Nubia. Our results indicate that despite the recent slab break‐off, Sinai has and is still being pulled (or pushed) toward the overriding Anatolia Microplate. The continued convergence possibly occurs because of a persistent slab pull force, a suction force induced by the down‐going detached slab and/or by the upper mantle flow induced by the Afar Plume.

Thermal Modeling of the Sanbagawa and Ryoke Belts

The Sanbagawa and Ryoke belts were formed in a convergent plate boundary along the eastern margin of Eurasia. Thermal modeling using the geological records of these belts as constraints allows quantitative estimates of both shear heating along the Wadati-Benioff zone and magma fluxes beneath the volcanic arc. In contrast to real-time observations of crustal movement and heat flow, rocks record changes in pressures and temperatures that occur over periods of several million years and can be used to examine conditions from the surface to the mantle. Thermal modeling combined with such geological records helps to bridge the gap in our knowledge between real-time observations of ongoing geological processes and the development of orogenies in convergent plate margins over geological time.

Rift obliquity in the Northern Volcanic Zone in Iceland using UAV-based structural data

Rift obliquity is known to be very frequent at divergent plate boundaries and shear components in the opening of fissures have been detected through seismicity in the monitoring of recent rifting events. However, the field structural evidence of obliquity in active rifts for Holocene events is still poorly studied, despite its relevance at both the lithospheric and upper-crust scales.We performed extensive high resolution UAV surveys on Holocene terrains in four areas along the active rift in the Northern Volcanic Zone of Iceland, to map structures and morphotectonic features in detail. We did a kinematic analysis and found an increase of the obliquity from N to S.We further discuss the origins of this obliquity and suggest that the combination of inherited structures reactivation, magmatic intrusions, and topography influence on the propagation of both fractures and magma intrusions can have an impact on the mismatch between rift opening direction and the far-field strain.We conclude that obliquity can be inferred for Holocene rifting events from surface structural data. It is thus fundamental to take into account the existing structural pattern and its (re)activation history in active rift zones to avoid missing out on possible preferential structure reactivation, magma propagation ways, and eruption sites during the monitoring of ongoing rifting events.

Path and Slip Dependent Behavior of Shallow Subduction Shear Zones During Fluid Overpressure

AbstractElevated pore fluid pressure is proposed to contribute to slow earthquakes along shallow subduction plate boundaries. However, the processes that create high fluid pressure, disequilibrium compaction and dehydration reactions, lead to different effective stress paths in fault rocks. These paths are predicted by granular mechanics frameworks to lead to different strengths and deformation modes, yet granular mechanics do not predict their effects on fault stability. To evaluate the role of fluid overpressure on shallow megathrust strength and slip behavior, we conducted triaxial shear experiments on chlorite and celadonite rich gouge layers. Experiments were conducted at constant temperature (130 and 100°C), confining pressure (130 and 140 MPa), and pore fluid pressures (between 10 and 120 MPa). Fluid overpressure due to disequilibrium compaction was simulated by increasing confining and pore fluid pressure synchronously without exceeding the target effective pressure, whereas overpressure due to dehydration reactions was simulated by first loading the sample to a target isotropic effective pressure and then increasing pore fluid pressure to reduce the effective pressure. We find that the effects of fluid pressure and stress path on the mechanical behavior of the chlorite and celadonite gouges can generally be described using the critical state soil mechanics (CSSM) framework. However, path effects are more pronounced and persist to greater displacements in chlorite because its microstructure is more influenced by stress path. Due to its effects on microstructure, the stress path also imparts greater control on the rate‐dependence of chlorite strength, which is not predicted by CSSM.

Linking 3D Long‐Term Slow‐Slip Cycle Models With Rupture Dynamics: The Nucleation of the 2014 Mw 7.3 Guerrero, Mexico Earthquake

AbstractSlow slip events (SSEs) have been observed in spatial and temporal proximity to megathrust earthquakes in various subduction zones, including the 2014 Mw 7.3 Guerrero, Mexico earthquake which was preceded by a Mw 7.6 SSE. However, the underlying physics connecting SSEs to earthquakes remains elusive. Here, we link 3D slow‐slip cycle models with dynamic rupture simulations across the geometrically complex flat‐slab Cocos plate boundary. Our physics‐based models reproduce key regional geodetic and teleseismic fault slip observations on timescales from decades to seconds. We find that accelerating SSE fronts transiently increase shear stress at the down‐dip end of the seismogenic zone, modulated by the complex geometry beneath the Guerrero segment. The shear stresses cast by the migrating fronts of the 2014 Mw 7.6 SSE are significantly larger than those during the three previous episodic SSEs that occurred along the same portion of the megathrust. We show that the SSE transient stresses are large enough to nucleate earthquake dynamic rupture and affect rupture dynamics. However, additional frictional asperities in the seismogenic part of the megathrust are required to explain the observed complexities in the coseismic energy release and static surface displacements of the Guerrero earthquake. We conclude that it is crucial to jointly analyze the long‐ and short‐term interactions and complexities of SSEs and megathrust earthquakes across several (a)seismic cycles accounting for megathrust geometry. Our study has important implications for identifying earthquake precursors and understanding the link between transient and sudden megathrust faulting processes.

Impacts of Groundwater Pumping on Subterranean Microbial Communities in a Deep Aquifer Associated with an Accretionary Prism

Accretionary prisms are composed mainly of ancient marine sediment scraped from the subducting oceanic plate at convergent plate boundaries. Anoxic groundwater is stored in deep aquifers associated with accretionary prisms and can be collected via deep wells. We investigated how such groundwater pumping affects the microbial community in a deep aquifer. Groundwater samples were collected from a deep well drilled down to 1500 m every six months (five times in total) after completion of deep well construction and the start of groundwater pumping. Next-generation sequencing and clone-library analyses of 16S rRNA genes were used to describe the subterranean microbial communities in the samples. The archaea: the prokaryote ratio in groundwater increased significantly from 1 to 7% (0 and 7 months after initiating groundwater pumping) to 59 to 72% (13, 19, and 26 months after initiating groundwater pumping), and dominant prokaryotes changed from fermentative bacteria to sulfate-reducing archaea. The optimal growth temperature of the sulfate-reducing archaea, estimated based on the guanine-plus-cytosine contents of their 16S rRNA genes, was 48–52 °C, which agreed well with the groundwater temperature at the deep-well outflow. Our results indicated that, in deep aquifers, groundwater pumping enhances groundwater flow, and the supply of sulfate-containing seawater activates the metabolism of thermophilic sulfate-reducing archaea.

Variable plate kinematics promotes changes in back-arc deformation regime along the north-eastern Eurasia plate boundary

The stretching of the lithosphere leading to back-arc basins formation generally develops behind arc-trench systems and is considered the consequence of slab retreat relative to the upper plate. Here, we examine the deformation regime evolution within the overriding plate due to subduction processes, using thermo-mechanical numerical simulations. We explore the north-eastern Eurasia plate boundary and the mechanisms of subducting Pacific plate since 57 Ma. During this time interval, several extensional basins formed along the Eurasia margin, such as the East China Sea, the Japan Sea, and the Kuril basin. Here, we increased the simulation complexity, with the inclusion of (i) the kinematic variability of the Pacific plate over the geological past with respect to a fixed Eurasia, incorporating time-dependent (i.e., temporally evolving) velocities computed from plate motion reconstructions; (ii) a Low-Velocity Zone within the asthenosphere, and (iii) a horizontal eastward mantle flow. Our results show a crucial role of the mantle flow for the development of lithospheric extension and back-arc basin opening, and a main kinematic control of the subduction trench position, which advances and retreats, into distance intervals in the order of ∼\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim$$\\end{document} 100 km, and providing stages of compression and extension in a back-arc basin.

Brittle faulting and tectonic stress history on the western margin of the Congo Basin between Kinshasa and Brazzaville: Implications for the evolution of the Malebo Pool and the Congo River

The history of brittle faulting and associated paleostress states has been investigated in the rapids of the Congo River at the outlet of the Malebo Pool between Kinshasa (Democratic Republic of the Congo) and Brazzaville (Republic of the Congo). The aim of this study was to unravel the Meso-Cenozoic tectonic evolution of the western margin of Central Africa, as well as to investigate the dynamic evolution of the Pool and the capture of the Congo Basin drainage system by the Lower Congo River through the Central African Atlantic swell that formed between the Congo Basin and the Atlantic Ocean. We reviewed the geological evolution of the Pool area using recent data and new field observations to develop a consistent stratigraphic model for both sides of the Congo River. We collected a significant amount of data faults and fractures in the early Paleozoic Inkisi arkosic sandstones, identified four successive data subsets and determined their related paleostress tensors trough an iterative process utilizing relative chronological indications. All paleostress states exhibit a a strike-slip faulting regime, with varying orientation of horizontal principal compression. The timing of faulting was constrained by indirect observations based on the updated understanding of the Pool's evolution. The geodynamic interpretation of the brittle stages considered potential sources of compressional stress arising from plate boundaries and from the oceanic ridge push. A polyphase brittle tectonic history was identified, characterized by the gradual appearance of newly formed faults and reactivation of the existing ones, leading to the gradual saturation of the host rock with fractures. The most recent event let to a significant imprint on the current landscape and is believed to have facilitated the connection of the drainage system of the Congo Basin to the Atlantic coast in the Oligocene period. The current stress state, determined from earthquake focal mechanisms, aligns with the most recent paleostress state identified determined from fault-slip data, displaying similar stress axis orientations but with consistent with a transpressional stress regime instead of a strike-slip one.

Temporary slip speed increases during short-term slow slip events with durations of one to three hours

AbstractTemporary slip speed increases with durations of 1–3 h were identified during short-term slow slip events in records of borehole and laser strainmeters in the Tokai region, Japan. They were found by searching for peaks of correlation coefficients between stacked strain data and ramp functions with rise times of 1 and 2 h. Although many of the strain steps were considered due to noise, some strain steps occurred with simultaneous activation of the deep tectonic tremors and shared source areas with the tremors. From 2016 to 2022, we observed five strain steps with simultaneous activation of tectonic tremors and coincidence of source locations with the tremors. Those strain steps occurred during short-term slow slip events and were temporary slip speed increases of the slow slip events. Those strain steps seemed to be related to successive occurrences with source migration of short-term slow slip events. The detrended strain steps corresponded to plate boundary slip events of moment magnitude around 5, which was consistent with the scaling law of slow earthquakes. Graphical Abstract

Intraplate stress distribution within the Indian Plate: Insights from finite element modelling

The Indian Plate is bounded by different types of plate boundaries like divergent, convergent/collisional, transform and consists of oceanic and continental lithosphere. Due to such tectonic settings, the intraplate tectonic stress distribution across the Indian Plate is very complex. Further, available in-situ stress information over the Indian Plate is relatively sparse; therefore, numerical modelling of maximum horizontal compressive stress (SHmax) is needed to comprehend regional stress patterns. In the present work, stress simulations are carried out to explore the effects of rigidity, lithospheric thickness, and boundary forces on the distribution of intraplate stress across the Indian Plate. The rigidity is derived from effective elastic plate thickness (Te), which is estimated through joint inversion of admittance and coherence using the wavelet transform approach. The stress field is simulated for four cases for different combinations of Te, lithospheric thickness and tectonic forces and compared with in-situ stress observations. The analysis of models shows that the magnitude of SHmax is primarily controlled by boundary forces, Te and scaled by lithospheric thickness; and its orientations are mostly governed by the applied tectonic forces and rigidity. Spatial distribution of earthquakes correlated well with the von-Mises stress and high von-Mises stresses are accumulated at earthquake prone regions. The stress field and von-Mises stress over major geological features are also analyzed and discussed. Apart from boundary forces, it is inferred that Te is one of the most significant factors which controls the stress field distribution.

Chronological, Petrogenetic, and Tectonic Significance of Paleoproterozoic Continental Crust within the North China Craton during the Global Tectono-Magmatic Lull

Abstract The 2.45–2.20 Ga period during the early Paleoproterozoic era is considered to have witnessed a global “Tectono-Magmatic Lull (TML)” and thus marks a relatively quiescent period. Our study unveils a 2.45–2.20 Ga magmatic suite from the Xiong’ershan area in the southern North China Craton, offering some key constraints on localized active tectonics during the TML. Zircon U-Pb dating shows Paleoproterozoic ages for the meta-basalt (2.31, 2.28 Ga), Na-rich meta-andesite (~2.33 Ga), tonalite-trondhjemite-granodiorite (TTG) gneisses (2.36, 2.30 Ga), K-rich granodiorite (~2.29 Ga), and monzogranite (2.33, 2.27 Ga). The meta-basalts geochemically and petrographically belong to calc-alkaline basalts and show distinctive Nb, Ta, and Ti contents and primitive mantle normalized patterns from different places in the Xiong’ershan area. Combined with their enriched εHf(t) values, the magmas were derived from subduction-related enriched mantle sources within a convergent plate boundary. The meta-andesites display high MgO content (average 4.5 wt%) and Mg# (44–57), strongly fractionated rare-earth pattern, calc-alkaline affinity, and negative Nb, Ta, and Ti anomalies. The TTG gneisses are of high SiO2 type (>62 wt%), high (La/Yb)N (17.5, 39.2), and Sr/Y (50.2, 104.3) and mostly display positive Eu anomalies and high-pressure type. Zircons from these rocks show a relatively narrow range of δ18O isotope values (5.35‰, 6.79‰) with εHf(t) isotope characteristics (−9.3, −3.3), suggesting derivation from partial melting of a thickened mafic lower crust. The youngest K-rich granodiorite and monzogranite show high K2O/Na2O ratios (0.65, 2.45). Variable molar ratio Al2O3/(CaO+Na2O+K2O) (A/CNK) and low zircon εHf(t) values suggest that the K-rich granitoids formed from the partial melting of different levels of crust. The presence of meta-basalt to andesite assemblages and diverse intermediate to felsic magmatic rocks implies magmatic activity within a convergent plate boundary tectonic environment with potential influence from plume-triggered extensional processes, supported by evidence of slab rollback and upwelling of mantle material.

Enrichment of arsenic in the Yarlung Tsangpo basin, Southern Tibetan Plateau: Provenance, process, and link with tectonic setting

Geogenic As enrichment is prevalent in freshwater systems worldwide, yet the connection between such enrichment and regional tectonics remains poorly understood, especially within orogenic environments like the Tibetan Plateau. To provide new insights into the sources and processes of arsenic, as well as its relationship with tectonics in the plateau basin, we selected the Yarlung Tsangpo basin, situated at the collision zone between the Indian (Himalaya block) and Eurasian plates (Lhasa block) as the study area. We analyzed and compared the chemical compositions of waters and sediments across various litho-tectonic domains. This also involved assessing arsenic concentrations in river waters, suspended particulate matter (SPM), sediments, geothermal waters, and rocks. The main findings are as follows: (a) High-arsenic geothermal waters and igneous rocks are common in the Lhasa block, and the discharge of these geothermal waters is the main source of dissolved As in river water, whereas erosion of igneous rocks contributes most of the As in sediments. (b) Metal oxides adsorb dissolved As, which leads to further As enrichment in the sediments. (c) The occurrence of high As concentrations in the southern Tibetan Plateau can be explained by a crustal cycling model. Arsenic likely originated in Neo-Tethys oceanic crustal rocks and was concentrated in the upper crust of the active continental margin (Lhasa block) during geological cycling at convergent plate boundaries. These studies highlight the role of tectonic settings on influencing arsenic enrichment in the plateau basin.

Discrimination between icequakes and earthquakes in southern Alaska: an exploration of waveform features using Random Forest algorithm

SUMMARY This study examines the feature space of seismic waveforms often used in machine learning applications for seismic event detection and classification problems. Our investigation centres on the southern Alaska region, where the seismic record captures diverse seismic activity, notably from the calving of marine-terminating glaciers and tectonic earthquakes along active plate boundaries. While the automated discrimination of earthquakes and glacier quakes is our nominal goal, this data set provides an outstanding opportunity to explore the general feature space of regional seismic phases. That objective has applicability beyond ice quakes and our geographic region of study. We make a noteworthy discovery that features rooted in the spectral content of seismic waveforms consistently outperform statistical and temporal features. Spectral features demonstrate robust performance, exhibiting resilience to class imbalance while being minimally impacted by factors such as epicentral distance and signal-to-noise ratio. We also conduct experiments on the transferability of the model and find that transferability primarily depends on the appearance of the waveforms. Finally, we analyse misclassified events and find examples that are identified incorrectly in the original regional catalogue.